Membrane tension lowering induced by protein activity 
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We present measurements of the fluctuation spectrum of giant vesicles containing bacteri- 
orhodopsin (BR) pumps using video-microscopy. When the pumps are activated, we observe a 
significant increase of the fiuctuations in the low wavevector region, which we interpret as due to a 
lowering of the effective tension of the membrane. 
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Membranes are self-assembled bilayers of surfactants 
or phospholipids, which form flexible surfaces. The me- 
chanical properties of membranes are essentially con- 
trolled by parameters, such as the membrane ten- 
sion, bending modulus and spontaneous curvature !l|. 
These parameters completely characterize the membrane 
fluctuation spectrum at thermal equilibrium, but they 
are not sufflcient to characterize real biological mem- 
branes, such as cell plasma membranes, which are non- 
equilibrium systems [2, ISj . Biological membranes are in 
general active in the sense that they are constantly main- 
tained out of equilibrium either by active proteins inside 
the membrane that are often ATP-consuming enzymes 
or by an energy flow, for example a lipid flux). Due to 
the complexity of active biological systems in vivo, recent 
studies have focused on in vitro biomimetic systems. An 
example of an active system is a giant unilamellar vesi- 
cle (GUV) which is rendered active by the inclusion of 
light-activated bacteriorhodopsin (BR) pumps. The bac- 
tcriorhodopsin pumps transfer protons unidirectionally 
across the membrane as a consequence of their conforma- 
tional changes, when they are excited by light of a spe- 
cific wavelength. Experiments on giant unilamellar vesi- 
cles containing bacteriorhodopsin pumps or Ca^+ATPase 
pumps have shown that the non-equilibrium forces aris- 
ing from ion pumps embedded in the membrane are 
able to significantly enhance the membrane fluctuations 
[J, |5|, la, |7| . In these experiments, the membrane tension 
has been fixed by means of micropipet aspiration and the 
corresponding excess area measured. The slope of such 
a curve defines an effective temperature, measuring the 
enhancement of the membrane fluctuations by the pro- 
tein activity. The amplification of the fluctuations due 
to activity has been originally predicted in Ref. Q . 

In this Letter, we report the first experimental mea- 
surements of an active fluctuation spectrum of GUVs 
containing bacteriorhodopsin pumps using video mi- 
croscopy. The details of the technique and of the analysis 
are given in Ref. [Qj. These measurements are comple- 
mentary to previous experiments using micropipets since 



they probe a range of wavevectors not accessible other- 
wise. Micropipet experiments provide information on the 
active fluctuations integrated over all wave vectors which 
require an independent calibration of the membrane ten- 
sion. On the contrary, video microscopy gives a direct 
measurement of the fluctuation spectrum from which the 
membrane tension can be extracted [9|. In this paper, we 
analyze theoretically the active fluctuation spectrum us- 
ing the model of Ref. [lO], which is an extension of the 
hydrodynamic model of Ref. @ to quasi-spherical vesi- 
cles taking into account the active noise of the pumps 
and localized active forces. 

Let us first discuss the passive behavior of giant vesi- 
cles. For the sake of simplicity we use notations appro- 
priate to quasi-planar membranes. We consider a quasi- 
planar membrane surface in the Monge gauge, in which 
it is defined by its height w(x) at position x. Two density 
fields are defined on the surface, i/' and tp K correspond- 
ing to the surface densities of proteins oriented with the 
two possible orientations [5|]. The membrane free energy 
is a function of the height field m(x) and of the imbalance 
of protein densities ipi^x.) = i'^ {'x) — ip-^{'x.). We expand 
the free energy to quadratic order in these variables: 

F[u,^Ij] = i f d^yi[K{V^uy + a{Vuf+xip^-2E'iljV^u], 

(1) 
where k is the bending modulus, a the surface tension, x 
a susceptibility coefficient and S a coefficient character- 
izing the coupling between the membrane curvature and 
the average orientation of the proteins. In the following, 
we neglect the curvature induced coupling; this is justi- 
fied for bacteriorhodopsin pumps since the incorporation 
of proteins does not lead to any measurable renormal- 
ization of the bending modulus [5|] . Note that for dilute 
proteins in the membrane, x = ksT /tiq, where n^ is the 
average density of proteins, of the order of lO^^m^^. In 
our experiments, only the vesicle contour at the equa- 
tor corresponding to a slice of the vesicle in the plane 
2/ = is recorded and is used to calculate the fluctuation 
spectrum: the measured quantity is 
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The passive equilibrium value of this fluctuation spec- 
trum is calculated using the equipartition theorem from 
Eas. mil 
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Fluctuation spectra are often analyzed in terms of power 
laws defined by gijix-^ cr, k) ~ q^. Two regimes can be dis- 
tinguished depending on the position of the wave vector 
qx compared to the cross-over wavevector qc = ^J a/ n. 
For q -^ qc, membrane tension dominates and v — —\, 
whereas for q ^ qc bending elasticity dominates and 
1/ = — 3. Our experiments lead to i/ ~ — 2, in the cross- 
over region between these two regimes. 

We now discuss active membrane fluctuations. In 
Ref. Q , a hydrodynamic theory was developed to calcu- 
late the non-equilibrium fluctuations of an active mem- 
brane containing ion pumps. This work has stimu- 
lated substantial theoretical interest, focused mainly on 
the general question of the proper description of non- 
equilibrium effects associated to protein conformation 
changes [ll|, [l^l • More recent developments of a similar 
hydrodynamic approach have lead to a general theoret- 
ical description of active gels [ij, [Ij, [T^ . In Ref. [lO] , 
the active force distribution is modeled as a superpo- 
sition of dipoles located along the membrane normal n 
(and not along z as in previous models), which generalizes 
Ref. Q . Experiments do not probe this force distribution 
directly but rather its first two moments. The first mo- 
ment represents the active contribution to the membrane 
tension and is denoted by Udip and the second moment 
Q represents the modification of the membrane bending 
moments by the activity of the force dipoles [iO| . Both 
active contributions to the tension and to the membrane 
curvature are also present in Ref. [l6| where the activity 
is due the the myosin molecular motors in the cortical 
layer bound to the cell membrane. 

The fluctuation spectrum of an active membrane in a 
quasi-spherical geometry can be calculated using an ex- 
pansion around a spherical shape with radius i?o to first 
order in the deviations from the sphere using spherical 
harmonics, with R ~ Rq + ni?o X); m '^imyim- The fiuc- 
tuation spectrum of a quasi-spherical vesicle is charac- 
terized by the amplitude of the spherical modes (|u;m|^) 
which are function of I only for symmetry reasons. In 
Ref. 0, the fluctuation spectrum of a quasi-spherical 
vesicle is compared to the corresponding spectrum for 
a planar membrane, and this comparison showed that 
for modes / > 5, the fluctuation spectrum of a quasi- 
spherical vesicle are very well approximated by that of 
a planar membrane at a wave vector g, (|M(g)p); at a 
wavevector q — 1/Rq, the two spectra are related by 
(|u(g)P) — i?Q(|u;,„p). For the sake of simphcity, we dis- 
cuss in the following the effect of the activity on the fluc- 
tuation spectrum only at the level of planar membranes. 
We have also simplifled the calculation of Ref. [l^, by 
neglecting the shot-noise, which is the intrinsic noise of 



the pumps [q| . This approximation is reasonable because 
the main effect of the shot-noise is to modify the fluctua- 
tion spectrum in the high wavevector region [i3| which is 
not the region where we observe a large modification of 
the spectrum due to activity in our experiments. We de- 
note here F2 the constant entering in Q = F21JJ [10]. This 
generalizes the notation of Ref. 5] , where this quantity 
represented the quadrupole moment of the force dipole. 
With these three approximations, we obtain the following 
spectrum 
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where D is the diffusion coefficient of the proteins which 
is of the order of lO^^^m^s^^ and 77 is the viscosity of 
water. Since the tension of the vesicles is at least of the 
order of 10~^N/m, we find that the term 4rjDq can be 
neglected with respect to a in the range of wavevectors 
relevant to the experiments [I3|. The effective tension 
a is defined as ct = tr -I- adip, where a is the passive 
contribution to the tension and <Tdip the active part due 
to BR proteins. After integration as in Eq. [3l we find 
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in which both terms on the right hand side contain ac- 
tive contributions when adip or F2 do not vanish. Note 
that Eq. |4] is very similar to the active spectrum given 
Ref. [5] with the correspondence F2 = 2w'Pa except for 
one important difference: the tension contains an active 
contribution but not the bending modulus whereas in 
Ref. Q it is the opposite. Because of this difference the 
active spectrum of Ref. [5| does not fit our experiments (it 
would predict a large effect of activity at high wavevec- 
tors whereas the effect is observed here mainly at small 
wavevectors). Another consequence of this difference, is 
that the active spectrum of Eq. [5] cannot in general be 
described in terms of an active temperature even in the 
low wavevector limit |l8l |. 

Experiments have been performed with GUVs made 
of Egg Phosphatidyl Choline at a molar fraction of 240 
lipids per protein, using the BR reconstitution protocol 
developed in Ref. [l0|. We have used a PEG passivated 
substrate and checked that vesicles fluctuate without ad- 
hering [9| . In Fig. m we show the measured fluctuation 
spectrum for the same GUV containing BR in active and 
passive states. To ensure reproducibility of the data, 
three consecutive passive and then, three active spectra 
were recorded. This flgure confirms that activity leads 
to an enhancement of the fluctuations as expected from 
previous experimental studies [5, 6]. It is important to 
point out that to observe this effect ImM of azide must 
be present in the solution. In the absence of this com- 
pound no enhancement of the fluctuation spectrum was 
observed [2^. This observation is consistent with sev- 
eral studies which suggest that azide enhances the pro- 
ton transfer in BR I2ll|. A large enhancement of the 
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FIG. 1: Fluctuation spectrum of a single GUV containing BR 
in a buffer with ImM sodium azide. The red (empty) symbols 
correspond to three consecutive recorded passive spectra, and 
the green (filled) ones to three consecutive active spectra. 
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FIG. 2: Fit of one active spectra of Fig. [T] using Eq. \5\ The 
parameter apparent k was fixed to 5.6 • 10~^®J obtained by 
fitting the passive spectrum. This fit gives a = 5.3 ■ 10~*N/m 
and Fz = 3.9 • 10"^** Jm. 



active as compared to passive fluctuations is observed at 
low wavevectors, which we attribute to a lowering of the 
membrane tension due to the activity. 

We first fit the passive spectrum using Eq. [3l This 
leads to cr = 3.9 • 10^^ ± 3 • 10~^N/m and an apparent 
bending modulus k = 5.6 • IQ-^^ ± 4.4 • IQ-^o J. The value 
of the bending modulus deduced from such a fit is over- 
estimated. Several artifacts perturb the measurements 
at high wavevectors: the pixel noise due to the discrete 
detection of the images, the integration time of the cam- 
era, and the effect of gravity Q. Taking into account 
these corrections at high q becomes very difficult in the 
active case, we thus have chosen a simple fitting pro- 
cedure, in which the corrections due to the integration 
time and to the gravity are neglected for both passive 
and active spectra. This procedure is justified by the ob- 
servation that both spectra are superimposed at high q. 
The determination of the bending modulus k is therefore 
delicate with our video microscopy technique while the 
measurement of the tension is reliable. With the value of 
the bending modulus kept fixed, we fit the three active 
spectra of Fig. [1] using Eq. [5l One fit which is shown in 
Fig. m leads to CT = 5.3 • lO^^N/m and Fa = 4 • IQ-^^Jm. 
Thus the active contribution to the tension adip ~ o — a 
is negative and of the order of —3.3 • 10~^N/m for this 
particular vesicle (see table I for a summary of the re- 
sults of the fits for the GUV of Fig. [1|). The lowering 
of the tension is systematic in all our experiments. The 
fact that (Jdip — — 3.3- lO^^N/m is very different from the 
estimate of 1.6 • 10^'^N/ni based on Refs. [a, 13 suggests 
that the modeling of the force distribution of these refer- 
ences is not appropriate to our experiments. In table I, 
we have also shown the exponent v, which is the apparent 
power law exponent of the spectrum. When the vesicle 
is active the exponent gets closer to —3, which is the 
expected value in the bending dominated regime. This 
is consistent with a lowering of the cross-over wavevec- 
tor qc — y/a/ K, and thus with a lowering of the tension 
since in this model k is not affected by activity. We have 
ignored any dependence of the protein activity on the lo- 



cal membrane curvature, an assumption justified by the 
observation that the main effect is a correction to the 
tension and not to the bending modulus. The value of 
the parameter F2 deduced from our fit is of the same or- 
der of magnitude as that estimated in Ref. @ based on 
micropipet experiments since the factor 2-3 increase in 
effective temperature measured in this reference leads to 
F2 ~ 9 • 10~^^ Jm. To summarize, Eq. 2] successfully de- 
scribes both the micropipet and the fluctuation spectrum 
measurements. 

We now discuss in more details the lowering of the 
membrane tension due to activity. We have checked that 
this effect is compatible with the constraints of constant 
surface and volume of the vesicle^ which must be im- 
posed if permeation is negligible |5|. To do so, we es- 
timate the excess area a — J q'^dq{\u{q)\'^)/4:/7T for pas- 
sive and active vesicles based on the above fluctuation 
spectra, by taking the lower bound of integration to be 
Qmin = 27r/i?o, with Rq — 10/im the radius of the vesi- 
cles; this lower bound corresponds roughly to the lowest 
measured value in Fig. [TJ In the accessible range of q val- 
ues, the excess area density q^{\u{q)\'^)/4:/TT of the active 
vesicles is always larger than the excess area density of 
the passive vesicles. At higher q wavevectors, it should 
be the opposite; but the q range where this would happen 
may be outside the accessible range. 

Another important question is the sign of the effect, 
which can not be flxed by symmetry arguments [4. In- 
deed, both orientations of the force dipoles @, [III, in- 
wards or outwards with respect to the membrane surface, 
are possible. The present experiments suggest that here 
the force dipoles must have an inward orientation (corre- 
sponding to contractile forces [14]) to produce a lowering 
of the tension. Another possibility would be to consider 
lateral dipoles. Here, we propose that the lowering of 
the tension is a consequence of electrostatic effects. This 
interpretation is not incompatible with the fact that BR 
conformational changes also play an important role. We 
have checked experimentally that a pH gradient builds 
up in our GUVs when they are activated, which suggests 



Spectrum 


a (x 10"'^ N/m) 


F2 (x lO"''** J.m) 
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passive 


3.9 ±0.3 


(imposed) 


-2.15 ±0.04 


active 1 
active 2 
active 3 


0.53 ±0.2 
0.79 ±0.06 
0.35 ±0.1 


3.9 ±1.1 
9.6 ±1.1 
7.2 ±0.6 


-2.7 ±0.04 
-2.78 ± 0.06 
-2.87 ±0.07 



TABLE I: Fitting parameters using the fluctuation spectrum 
of Eq. [5] The passive spectrum is fitted with the condition 
F2 = 0. Then three consecutive active spectra are fitted with 



the same vesicle using the value of Ka 



5.6 • 10" 



±4.4 



10~ J determined from the fit of the passive spectrum. From 
this, the membrane tension ct, F^ and the exponent v are 
deduced. 



that a significant voltage drop exists across the mem- 
brane. Furthermore in membranes containing BR, the 
activity of the BR can be suppressed by applying an ex- 
ternal voltage drop across the membrane [22| . In view 
of these observations, we propose that Maxwell stresses 
created by ion transport lead to a renormalization of the 
tension as shown in a simple model developed by two of us 
[23] . In this model, the electrostatic correction to the ten- 
sion is proportional to the square of the electric current 
going across the membrane. Maxwell stresses are large 
because of the large mismatch in dielectric constants at 
the membrane, and they tend to reduce the membrane 
area, thus producing a negative electrostatic contribution 
to the tension. A lowering of the membrane tension due 



to the application of normal electric fields very generally 
leads to instabilities of the membrane [2J, [2^ . 

To summarize, we have presented what are to our 
knowledge the first measurements of active fluctuation 
spectra of GUVs containing BR proteins using video 
microscopy. The experiments show that activity en- 
hances the fluctuations. The effect is strong in the low 
wavevector region, which we interpret as a lowering of 
the tension due to activity. This effect could not have 
been detected in micropipet experiments which are done 
at constant membrane tension. A possible candidate for 
explaining the lowering of the tension would be that it is 
caused by Maxwell stresses due to the transport of ions 
across the membrane. Further experiments (using for 
instance patch-clamp techniques with the same GUV or 
other GUVs containing ion channels rather pumps) are 
required to confirm these suggestions. 
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